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Abstract. Generic wave dislocations (phase singularities, optical vortices) in three 
dimensions have anisotropic local structure, which is analysed, with emphasis on the 
O [ twist of surfaces of equal phase along the singular line, and the rotation of the local 

anisotropy ellipse (twirl) . Various measures of twist and twirl are compared in specific 
examples, and a theorem is found relating the (quantised) topological twist and twirl 
for a closed dislocation loop with the anisotropy C line index threading the loop. 



PACS numbers: 42.25.Hz, 03.65.Vf, 02.40.Hw 
1. Introduction 
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Q\ • Phase singularity lines in 3-dimensional scalar waves, that is, nodal lines or wave 

\o 

O 

G^ ■ or energy, flow. Being places where the phase is singular, all phases are found in the 
^ . vortex neighbourhood, and the phase change is arbitrarily fast there. The surfaces on 
O . which phase is constant all meet on the singularity line, and are often twisted around 
^ ■ the singularity as helicoids, as with the familiar screw dislocations PP; it is this twisting 
* ^ . that distinguishes screw dislocations from edge dislocations, whose phase structure does 

not rotate along the dislocation line. 
CLc In two dimensions, phase singularities are generically points where the intensity 

of the wave vanishes and the phase is singular. Knowledge of the phase contour lines 
^ . near the point singularities can be sufficient to fill in the phase structure of the rest 
. of the field. In 3-dimensional fields, topological singularities form a network of lines in 
space, often described as a 'skeleton' of the spatial pattern |21 EI - in complex scalar 
fields, the phase structure local to the phase singularity lines show how the global 
phase field is constructed upon this skeleton. The simplest 3-dimensional property 
of a phase singularity is its sense, or topological current: the direction endowed on 
the singularity line by the right-handed increase of phase. The topological current 
direction is preserved along the singularity line. The sense of the phase helicoids twisting 
around the singularity is independent of this; for instance (as demonstrated below), if 
a dislocation is embedded in a plane wave whose propagation direction is parallel to 
the topological current direction, the phase surfaces ending on the singularity form left- 
handed helicoids; if the propagation is antiparallel to the dislocation sense, they are 
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(a) 



(b) 




Figure 1. Two surfaces of constant phase (mod tt) in the vicinity of a twisted wave 
dislocation, (a) Uniform twist. This is a surface of constant phase (mod tt) of the simple 
screw dislocation (|1.2ll with k — (0,0, fc), and is a uniform helicoid. (b) Nonuniform 
twist. This constant phase surface is an example of the generic situation, calculated 
from with /3 = 2a and c = -1/2. 



right-handed. The twist is rather more comphcated in more general situations, such as: 
when the singularities evolve off-axis in gaussian beams, with a noncanonical transverse 
shape in EI; when they are knotted, linked, or braided [HIIZIIH]; in isotropic random plane 
waves, when the singularities are tangled in a nontrivial way j^]; when they form the 
characteristic 'antelope horns' of the elliptic umbilic diffraction catastrophe ^Oj- The 
zero contour surfaces of the real or imaginary parts of the field are often used to locate 
phase singularities and describe their geometry OIH]; these are, of course, special cases 
of the general phase contours. 

My aim here is to describe this 3-dimensional phase twist structure in the general 
case, illustrated by simple examples of dislocated waves. As these examples show, it is 
very easy to find dislocations in waves that do not have a uniform phase screw structure, 
and this situation is, in fact, what one would expect generically. The main complication 
to the description comes from the fact that different phase helicoids twist at different 
rates, and this must somehow be averaged to give an overall sense of dislocation twist 
at a point on the singularity line. This is related to the fact that the local phase 
structure transverse to a dislocation line is generically anisotropic, squeezed into an 
ellipse in El El 1121 Uni mi, complicating the averaging around the dislocation. An 
example of two phase helicoids near a dislocation are shown in figure ^ the first has a 
uniform helicoidal structure, the other does not. 

The notation of [HIEI will be followed. In particular, = ip{r) denotes a complex 
scalar wave, and in terms of real and imaginary parts, or amplitude and phase. 



Time dependence will not be considered here, and terms such as twist and twirl refer to 
variation with distance, not time. Local cartesian coordinates r = {x, y, z) or cylindrical 
coordinates i?, 0, z will be used, with the dislocation along x = ?/ = r = 0, and phase 
X increasing with (that is, the topological current is parallel to the +2;-direction). 
Important formulae will be given in coordinate-free form, with the dislocation tangent 



^ = ^ + i77 = pexp(ix). 
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denoted by T (which is (0,0,1) in local coordinates), and the directional derivative 
along the dislocation by (which is dz* in local coordinates). Several of the results 
described here were also discussed in jT2], particularly in section 2.8, pages 54-58. 
The epitome of a twisted wave dislocation is in the example |lj 



This wave is a local solution of the Helmholtz equation V^V + ^^"^ = 0- It is a screw 
dislocation when k = (0,0, k), and each of the surfaces of constant phase is a regular 
helicoid with pitch 27r//c, reminiscent of the atomic planes in a crystal lattice near a 
crystal screw dislocation jT^I. It may be compared with an edge dislocation, where in 
()1.2|1 . k = (0, fc,0). In this case, the phase structure does not change with z, and the 
phase lines have a structure similar to the arrangement of the atomic planes near a 
crystal edge dislocation A general k in ()1.2j) gives a mixed edge-screw dislocation. 
Far from the dislocation, the phase of ()1.2|1 has a plane wave character with wavevector 



This example, and its crystal analogy, motivated Nye [T^ to define the Burgers 
vector of a wave dislocation 



where subscripts after brackets denote quantities held constant, other subscripts denote 
partial derivatives, and coordinates are local. The dislocation ()1.2|1 has Burgers vector 
b = k, which is the wavevector of the asymptotic plane wave. This gives the expected 
results of the vector being parallel and perpendicular to the screw and edge dislocation 
lines respectively. However, there is a danger of taking the analogy between wave 
dislocations and crystal dislocations too far, and there are problems with this definition 
of a Burgers vector for a wave dislocation in other rather general situations. 

In the derivation of ()1.3p in ^^1; it is assumed that the dislocation is curved and 
moving (i.e. the wave is narrow-band); for monochromatic waves, the dislocation is 
stationary, and many dislocation lines occurring in optics are straight lines (often parallel 
to the beam direction). If the only available information about the field is local to the 
dislocation, for instance, a finite Taylor expansion of ()1.2|1 about x = y = 0, bz is not 
defined. In more general cases of wave interference, such as isotropic random plane wave 
superpositions [Hj, there is no overall propagation direction, and only local properties of 
the dislocation, determined using derivatives of ip on the dislocation, are relevant. 

It is desirable to have a measure of screwness, determining the twist of the local 
phase structure along the dislocation, in terms only of local derivatives. The following is 
an exploration of this twist geometry, and its realisation in simple beam superpositions. 

2. The geometry of twist and twirl 

The twist Tw is defined to be the rate of rotation of phase along a dislocation, according 
to several different measures to be described. The simplest geometric twist is the rotation 
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of a surface of constant phase along the dislocation, which is also the rate of rotation 
of the normal to the phase surface. For convenience, a surface of constant phase in the 
vicinity of the phase singularity will be called a phase ribbon [TTf I18j: for most of the 
discussion, only the ribbon geometry is important. Each plot (a) and (b) of figure 
shows two phase ribbons, with phases differing by vr. 

Fixing a particular phase Xo? the normal to the Xo-ribbon is 

Uxo = Re{V?Aexp(-ixo)} 

= V^cosxo + Vr/sinxo- (2.1) 
Using as a local azimuthal coordinate, the rate of twist Tw{xo) of the Xo-ribbon is 

Twixo) = i(l>z)x=xo 

= dz arctan ^"'^ 
U 

- (2.2) 

XO 

Applying this formula to the wave ()1.2p gives —kz. Apart from the sign, this is the 
z-component of the Burgers vector. The — sign originates in the fact that the helicoid 
is left-handed, a general property of dislocations whose topological current is parallel to 
the propagation direction [T21 E] • 

The phase structure transverse to the dislocation in (jl.2p is isotropic, and all the 
phase ribbons twist at the same rate as each other, and along the dislocation (one such 
is shown in figure Q (a)). General dislocations are anisotropic, however, and the local 
contours of intensity are elliptical, with a corresponding squeezing of phase gradient 
Vx- It can be shown 011211111 that {x<t>)z=Q = R^oj/ p^, where u is the vorticity on the 
dislocation, defined as u = |V^ x V?7| = iVip* x Vip\/2 (the direction of this vector 
gives the topological current). The plots in figure El indicate this aspect of elliptic phase 
squeezing. In fact, both the ellipse and the phase squeezing are accounted for by the 
gradient vector Vip on the dislocation. As with all complex vectors, it is associated 
with an eUipse, traced out by as Xo changes; this ellipse is the same shape as that 
described by p^ and Vx, but with axes exchanged. 

It is more natural to quantify the twist of the entire singularity core, rather than 
simply for a fixed phase. Therefore, it is necessary to average Tw{xo) over all phases, 
although it is not clear whether averaging with respect to the phase x or azimuth is 
appropriate. It is possible to calculate both averages: the phase average (integrating 
Tw{xo) with respect to xo) has a complicated form involving the axes of the phase 
ellipse, and is in jl2j equation (2.8.6); the azimuth average (integrating Tw(xo) with 
respect to (p) has a simpler form, and is ([121, equation (2.8.7)) 

T.^^f rd,r»w^M^^^:2iW}. ,,3, 

Z-n Jq Zuj 
This depends only on derivatives up to the second of on the dislocation line. As 
expected, it gives —k-^ for p.2|) . 
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Berry jinj chose to average phase instead by examining the rate of change of phase 
at a fixed azimuth {xz)<j>=4,o- The result of a particular averaging, he defined the screwness 
a to be 

-Im{VV;* ■ VV;^ 

^ = \vW ■ ^ ^ ^ 

(o" defined here is the negative of that defined in [IH]-) The screwness for fll.2|) is —kz- 
The difficulty in defining the total twist arises because the different phase ribbons 
twist at different rates, due to the phase anisotropy ellipse associated with Vip. Along 
the dislocation line, the anisotropy ellipse itself may rotate, as well as change its size 
and eccentricity. Because the rotation of the ellipse is independent of the phase twist, 
it will be referred to as the twirl tw of the dislocation line, and may be found as follows. 

The complex vector field Vip shares geometric features associated with vector 
polarization fields [T^ IT^ ^. In particular, in local coordinates, V^^ is confined to 
the xy-plane, and therefore parameters, describing all the geometric properties of the 
ellipse, may be defined: 

S2 = my + s, = -iimy + m;) = ^u. (2.5) 

These parameters describing the anisotropy are analogous to the Stokes parameters 
in polarization; the anisotropy ellipse is related to and these parameters as the 
polarization ellipse is related to the electric vector and the Stokes parameters. The 
parameters ()2.5|1 do not, themselves, have anything to do with polarization. The 
anisotropy ellipse, depending only on Vip, may be defined at any point of the scalar 
field, not only on a dislocation; the ellipse is circular (or linear) generically along lines 
in space [HI 123 1121 • 

The azimuthal angle of orientation of the major ellipse axis is arg(S'i + 182)/ 2. 
The twirl tw^ may therefore be defined as the rate of change of this angle along the 
dislocation line: 

1 SiS'2 - S2S[ 

The denominator is equal to Sq — 5*1 = |V^/'|'^ — Auj'^, which is zero when the ellipse is 
circular and twirl is not defined. 

The natural measure of phase around the ellipse associated with a complex vector 
is the rectifying phase Xr [121 1211 (^^^ polarization ellipses, it also called phase of the 
vibration [2]). Xt is defined such that the complex vector exp(— iXr)V^ has orthogonal 
real and imaginary parts (the real part along the ellipse major axis, the imaginary along 
the minor), and can be shown [21J to be equal to aigCVip ■ Vip)/2. The phase twirl tw^ 
may be defined as the rate of change of this phase along the dislocation: 

tz., = -(^arctan^^^-^J . (2.7) 

I The two fields do not have identical structures, since free field vector solutions of Maxwell's equations 
are divergence free, whereas ^ip, being a gradient field, is curl free. 
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Figure 2. The sequence of transverse phase hnes (separated by an equal phase 
difference 7r/6) and anisotropy elhpse of a twisting, twirhng dislocation of the form 
(123, with /3 = 2a,c = -f/2. the 2;-spacing between each frame is -k jla. The twist 
and twirl have opposite senses. The surface in figure (b) is swept out by a pair of 
opposite phase lines in this figure; the parameters in the two figures are the same. 



This gives a natural measure of the rate of change of phase with respect to the elhpse 
axes. Its form is not particularly simple when the derivative in fl2.7j] is taken, although 
it is easily seen that the denominator is | V"?/^!^ — 4a;^. This implies that tw^ is not defined 
when the ellipse is circular (isotropic). 

The two twirls here defined in ()2.(ij) . ()2.7|1 may be combined to give a new measure 
of the phase twist. Since tw-^ measures the rate of change of phase with respect to the 
ellipse axes, its negative gives a sense for the helicoid phase twist with respect to the 
ellipse. Therefore, the difference tw^ — tw^ gives an ellipse-defined twist Twew, which 
can be shown to be 

Re{T-(VV;*AVV;0+iVV;'--VV;^} 
T^eii = tw, - tw, = |vv;p + 2c. • ^^-^^ 

Although neither type of twirl is defined when the ellipse is circular, Twew is, and for 
the example ()1.2|) . it is —k^, as desired. It is also interesting to note that Tweii is the 
sum of the numerators of Tw^ and a, divided by the sum of the denominators. 

A simple example of a twisting, twirling dislocation is found in the sum of two screw 
dislocated waves, with dislocations in opposite directions, and with different pitches 
2n/a,2n/(3 : 

ipi = {x + iy) exp(ia;z) + c{x — iy) exp(i/3z) (2.9) 

with c in general complex, with 1 > |cp ensuring that the dislocation along x = y = Q 
has topological current in the +2;-direction. 
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(|2.9|) could represent, for example, the sum of two copropagating but 
counterrotating order one Laguerre- Gauss or Bessel beams with different components, 
in the vicinity of the 2;-axis. A possible realisation, in terms of Bessel beams, is 
exp(i0 + iaz)Ji{Vl-a^R) + cexp(-i0 + i/3^) Ji(^l - (3^)^(1 - a^)/{l - (3^). The 
rotation of a tranverse interference pattern along a beam, achieved though superposing 
beams different axial phase dependencies, was also studied by Courtial where the 
Gouy phases of two superposed singular beams were different; here, the components 
of the two dislocated fields are different. 

The phase change around the dislocation in ()2.9p is anisotropic for c 7^ 0; the 
Stokes parameters Sq, S3 are constant, indicating the anisotropy ellipse has a constant 
area 7r(l — |cp) and eccentricity 2|c|^/^/(l + |c|). ()2.9|) is therefore a local normal form 
for dislocations whose twist and twirl are much greater than their curvature or rate of 
change of anisotropy. 

The two types of twirl are computed to be 

tw^ = -[a- (3)/2, tw^ = {a + (3)/ 2 for (2.10) 
The various rates of twist are 

Tweii = -a I 

Tw^ = -{a - |c|2/?)/(l - IcH \ for (2.11) 
a =_(a + |c|2/5)/(l + |cp) J 

For this example of a dislocation with a uniform twirl, and constant anisotropy, 
the different measures of twist are, in general, different. They are equal if c — > 0, 
(yielding the screw dislocation in ()1.2j) ). or if the twirl is zero, i.e. a = (3 (this was the 
case discussed in |ll2j). Of the different measures of twist, it appears that Twew agrees 
numerically the most with intuition: if /3 is a positive integer multiple ma of a, the 
ellipse will undergo m — 1 rotations as the phase undergoes one, and the entire pattern 
is periodic (as in figure|21), with period 2n/a. Tweii is the only measure of twist to reflect 
this. Equivalently, it can be argued that since the first term on the right hand side of 
()2.9|) defines the dislocation direction, it also defines the twist, and the effect of the 
second term is merely to modulate the pattern to produce the twirl. 

3. Closed dislocation topology 

When phase singularity lines form closed loops, certain topological identities must 
be satisfied. In relation to twist, continuity of the wavefunction requires that the 
total number of twists of each phase ribbon must be a (positive or negative) integer, 
which is the same for each ribbon. This screw number is therefore a property of the 
dislocation loop, and is positive if the topological twisting is right handed with respect 
to the dislocation direction, negative if left handed. For obvious geometric reasons, a 
dislocation loop with nonzero screw number will be called a closed screw dislocation. 

The importance of the screw number is that it gives the dislocation strength 
threading the loop, by the twisted loop theorem: the screw number m of a strength 



Local phase structure of wave dislocation lines 



8 




Figure 3. A phase surface (mod tt) of the wave V'cioscd with m — 1. A straight vertical 
dislocation line threads a closed screw dislocation, and the local phase ribbons join up 
in the form of a noncompact torus. 

1 dislocation loop is equal to minus the dislocation strength threading the loop (in a 
right handed sense). This result is discussed and proved in [T21 CHI IIHI 1211 IE]- If the 
dislocation loop is planar, then the integral of the ribbon twist around the loop divided 
by 2tt gives the screw number. If the loop is nonplanar, then the Calugareanu- White- 
Fuller theorem |23| implies that the writhe of the curve must be added to the twist 
integral. Only planar curves will be considered here. 

A simple wave containing a closed screw dislocation can be made from a 
combination of polynomial waves in cylindrical coordinates (7j: 

^closed = i?'""' exp(im0) exp{ikz){R^ - Rl + 2i(|m| + l)z/k). (3.1) 

This wave has a closed screw dislocation in the z = plane at R = Rq, with screw 
number m, its topological current directed opposite to the increase of (p. It is threaded 
by a strength m dislocation up the z-axis (its sense in +z). In jHl Ej, high-strength 
loops with similar geometry were found. If m = 0, the loop is the familiar closed edge 
dislocation loop. Tw^,Twe\\ and a are —m/Ro on the closed loop; the screw number is 
m. The twirl is zero for this loop, and the anisotropy ellipse axes are oriented in the 
R, z directions. 

The nature of the twisted phase ribbons near the dislocation have consequences for 
the global topology of the total phase surface, that is, the wavefront. Figure El shows 
a surface of constant phase (modulo vr) of the wave V'ciosed with m = 1. This surface is 
a 'noncompact torus' (just as the plane is a 'noncompact sphere'), extended to infinity 
because of the infinite straight dislocation on the axis. Unlike a compact torus, there is 
no way of distinguishing the two sides of the surface. The discussions in [T21 UHl show 
that such complex wavefronts are inevitable with closed screw dislocations. 

The closed dislocation of (jH.lll has zero twirl. Waves with closed twirling 
dislocations loops can be made by superposing waves of the form '?/'ciosed with dislocations 
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of opposite sense and different topological twists: 

^2 = exp(imi0) exp(ifcz)(/?2 - Rl + 2i(|mi| + l)z/k) 

+ cE}'^^\ exp(im20) exp(-iA;z)(i?2 -B^- 2i(|m2| + l)z/k). (3.2) 

When mi,m2 are different, we may expect the twirl to be nonzero, as with ()2.9|1 . 
Not every combination of mi,m2 yields the expected twirling wave however, because 
the complicated threading interference structure is affected in a nonlinear way. For 
simplicity, k and Rq are taken as 1. 

Choosing mi = 1, m2 = 0, c = 2/3 gives a wave with a twirling dislocation, and the 
measures of twist are 

TW ell = -1 1 

Tw^ = -3(6 + cos0)/(2(7 + 3cos0)) > for ^/'2 with mi = 1, m2 = 0, c = 2/3.(3.3) 
a = -3(15-2cos0)/(53-12cos0) J 

The screw number of this wave is —1, since the phase structure of the summand with 
larger coefficient dominates. Only Twdi integrates to this integer; the others give 
irrational numbers, and this is the case for other choices of mi,m2,c. The integrability 
of Tweii may be explained by the fact that it is the derivative of the difference of angles 
()2.8|) . the others are averages of an angle, and taking the average does not commute 
with integrating around the dislocation loop. Thus, of all the various twists considered, 
only Tweii can be used as a topological twist. 
The twirl is 

t^/;^ = -(8 - 6cos0)/(25 - 24COS0) for ^2 with mi = 1, m2 = 0, c = 2/3, (3.4) 

and, by ()3.3|) and ()2.8|) . tw^ = 1 + tw^. The topological twirl around a closed loop is 
also quantised, equal to the number of rotations of the anisotropy ellipse around the 
loop, although it only needs to undergo a half turn to return to itself smoothly p| I21j. 
Integrating tw^ around the loop gives a topological twirl of —1/2 : with respect to 
the anisotropy ellipse, the loop is a Mobius band. The sense of topological twirl and 
topological twist are independent. 

A natural question to ask is whether the topology of twirl gives an analogue of the 
twisted loop theorem. In fact, it does: the ellipse rotation around a loop is related to the 
anisotropy C line index enclosed by the loop. Anisotropy C lines, where the anisotropy 
ellipse is circular, occur when Vip ■ is zero; they correspond to phase singularities of 
this complex scalar, whose phase is twice the rectifying phase Xr- The anisotropy C line 
index therefore is half this phase singularity strength, and is equal to the integral of tw^ 
around the loop, divided by 27r. This fact, together with the twisted loop theorem and 
the definition ()2.8|) of Twew-, leads to the twirling loop theorem, which may be stated as 
follows. 

The anisotropy C line index threading a closed dislocation loop (in units of 1/2) is 
equal to the topological twirl minus the topological twist of the loop, i.e. the number of 
rotations of the anisotropy ellipse around the dislocation, minus the number of rotations 
of the phase structure, in a right handed sense with respect to the dislocation strength. 
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Figure 4. The probability density functions of the different twist and twirl measures. 
In ascending order of y-axis interception, they are: twirl tw (thick black line) , helicoid 
twist Tw{xo) (thin black line), azimuth-averaged twist Tw^ (dotted line), screwness a 
(dashed line) , and ellipse-averaged twist Twcw (thick grey line) . The distributions are 
all in units of characteristic twist Twc- 

This gives a topological role to anisotropy C lines. Anisotropy L lines, where uj = 
(that is, where all the phase surfaces share a common normal), govern the reconnection 
of dislocation hnes [71 123 1211 • 

4. Twist and twirl in isotropic random waves 

As a final example, twist and twirl are considered in 3-dimensional isotropic random 
waves, that is, superpositions of plane waves with isotropic random directions and 
phases, whose ellipse anisotropy statistics have been calculated [111121 • These random 
waves might occur in monochromatic waves in a large, chaotic cavity, or (a scalar 
caricature of) black-body radiation; the wave dislocations form a complicated tangle. 
Nevertheless, this random wave model has many statistical symmetries, the averages 
being spatially and temporally invariant. 

The details of the calculations are omitted, but the calculations are reasonably 
straightforward, following the methods of [01 Ej, in which statistics of many geometrical 
properties of dislocations were calculated. The calculations proceed by taking advantage 
of the fact that ip and its derivatives have gaussian distributions, the details depending 
on the power spectrum of the waves considered. For all of the probability density 
functions, the twists and twirls are in units of the characteristic twistTwc = {kA/bk2Y^'^ , 
where fc„ is the nth moment of the power spectrum; for monochromatic waves of 
wavelength A, this is 27r/v^A. 

All measures of twist and twirl discussed here {Tw{xo) in fl^-^P - Tw^ in (j2.3|) . a 
in ()2.4|) . tw^ in ()2.6p . tw^ in ()2.7p . and Tw^w in ()2.8p ) involve some combination of 
first derivatives of the dislocation with their dislocation directional derivative, divided 
by another function involving first derivatives only. In the calculation, the gaussian- 
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distributed second derivatives are integrated first, giving a function involving uj and 
G = iV'j/'p, wliicli is tlien integrated using tfie probability distributions described in 

PIE]. 

The probability distribution functions of the various twists and twirls are found to 

be 

3 1 



4(1 + ^2)5/2' 

35 

32(l + t2)9/2^ 



Tw. 



'ell 



1 + t2)7/2 J 

V2(ll + 8t2 + 32t4 



;i + 2t2)7/2 (l+t2)3/2 



Ptu,{t) =-— — — ^ (2 + t2)(3 + 3t' + 8t^)E 



32t4(l + t2)5/2 ' ' ' \l+t 

2{^ + ^e + 2t')K(-^\], (4.1) 



where K represent the complete elliptic integrals of first and second kinds These 



are plotted in EJ they are all symmetric, since right-handed and left-handed screw 
dislocations are equally weighted in the ensemble. All of the twist distributions have 
power law tails. The fluctuations of the Ptw are the largest (the second moment does not 
converge). The fluctuations of Twew are the smallest (second moment is logv^ — 1/4), 
providing further support for the preference of Twew- 



5. Discussion 



Twist and twirl for anisotropic 3-dimensional wave dislocations are not important for 
the simplest optical vortices because they are isotropic and therefore not subject to 
the subtleties described here. However, it is important to note that the screw-edge 
distinction only applies in special cases such as (jl.2|) . and in more chaotic fields, only 
local phase geometry may be appealed to. Various measures of twist and twirl were 
introduced in |2l their topology was examined in section |21 and statistics in HJ The 
conclusion from each of these deliberations was that the ellipse-defined twist Tw^w fl2.8p 
has the most desirable properties of the various measures. 

Experimental verification of the twirling loop dislocation (j2.9p should be 
straightforward; with measurements of the phase, the twist and twirl structures of more 
complicated, curved dislocations may be found, although the full experimental analysis 
of the twist and twirl of an arbitrary dislocation line is likely to be rather difficult. 

It is natural to ask whether there is any curvature structure (i.e. second spatial 
derivative of ip) transverse to a dislocation, providing a transverse analogue to twist. 
This would provide the geometric counterpart to the transverse component of the 
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Burgers vector, as twist is related to the longitudinal part. For the standard edge 
dislocation ()1.2|) . the Burgers vector is normal to the phase surface (mod it) whose local 
transverse curvature is zero (a phase saddle in the transverse plane also occurs on this 
surface); it is tempting to guess this may be the desired structure. After some analysis, 
however, it can be shown that there is no unique phase contour whose local transverse 
curvature vanishes: generically, there are either one or three phase contours with this 
property. The number of phase contours in question is given by the number of real roots 
of a certain cubic in tan 0, reminiscent of the number of straight lined terminating on a 
line field singularity (lemon versus star and monstar) |2Hll^- An example of a solution 
to the wave equation with three phase lines of vanishing transverse curvature is 



Phase patterns for dislocations in two dimensions with this property are plotted in |29j . 
figure 6 (f). It is not certain, therefore, whether local dislocation geometry can give a 
transverse direction to a Burgers vector. 

The philosophy of studying phase singularities is that they give information about 
the global phase structure: they are a 3-dimensional 'skeleton' for the entire field. Twist 
is an important, intrinsic property of phase singularities, and shows how this information 
may be gleaned from the local dislocation morphology; closed screw dislocation loops 
indicate nontrivial wavefront topology, as in figure El Phase anisotropy gives rise to 
further structures, such as twirl, and anisotropy C and L lines. A possibility of a further 
extension is analysing the gradient of the anisotropy scalar Vijj ■ Vip itself: its phase 
singularities (the anisotropy C lines) themselves can be twisted, twirling, and therefore 
related to even higher singular morphologies. The twirling loop theorem shows how the 
first members of this hierarchy of anisotropy structures are coupled. 
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